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ABSTRACT 

We investigate the relationship between the present day optical luminosity 
function of galaxies and the X-ray luminosity function of Seyfert Is to determine 
the fraction of galaxies which host Seyfert 1 nuclei and their Eddington ratios. 
The local type 1 AGN X-ray luminosity function is well reproduced if 1% of 
all galaxies are type 1 Seyferts which have Eddington ratios of ~ 10"'^. However, 
in such a model the X-ray luminosity function is completely dominated by AGN 
in E and SO galaxies, contrary to the observed mix of Seyfert host galaxies. To 
obtain a plausible mix of AGN host galaxy morphologies requires that the most 
massive black holes in E and SO galaxies accrete with lower Eddington ratios, or 
have a lower incidence of Seyfert activity, than the central black holes of later type 
galaxies. 

Key words: galaxies:active - galaxies: luminosity function - accretion 



1 INTRODUCTION 

Accretion onto a massive black hole has remained the stan- 
dard paradigm for active galactic nuclei (AGN) for sev- 
eral decades (eg Lynden-Bell 1969, Rees 1984). In recent 
years, a considerable amount of evidence has accumulated 
for the existence of massive black holes in many, perhaps 
all, galaxies (Magorrian et al. 1998, Kormendy & Rich- 
stone 1995). The connection between the formation and 
evolution of massive black holes and their host galaxies 
is currently the subject of great interest and investigation 
(eg Fabian 1999, Salucci et al. 1999, Cattaneo, Haehnelt 
and Rees 1999). 

The luminosity function of AGN has been used for 
many years to track the statistical evolution of the AGN 
population with cosmic epoch, but does not allow distinc- 
tion between a relatively small population of long lived 
AGN or many short lived generations of AGN. The Ed- 
dington ratio, i.e. the ratio of an object's luminosity to 
its Eddington luminosity, was proposed to be a power- 
ful discriminator of AGN activity patterns by Cavaliere 
& Padovani (1988) with long lived AGN having low Ed- 
dington ratios (~ 10"*) and short lived AGN having high 
Eddington ratios (~ 1). In a variety of wavebands, from 
radio to X-ray, the luminosity function of AGN has a char- 
acteristic two-power law shape with a knee dividing the 
low and high luminosity objects (eg Dunlop & Peacock 
1992, Boyle, Shanks & Peterson 1988, Page et al. 1997). 
It has an obvious resemblance to the luminosity function 
of galaxies, which is not surprising since AGN are found in 



galaxies. Only a few attempts have been made to relate the 
galaxy and AGN luminosity functions, predominantly in 
the form of models for the joint formation of galaxies and 
AGN (eg Kaulfmann & Haehnelt 2000, Monaco, Salucci, 
& Danese 2000 and Haehnelt & Rees 1993) yet the direct 
comparison of galaxy and AGN luminosity functions can 
shed light on two fundamental properties of the black hole 
population: the fraction of massive black holes which are 
active, and their Eddington ratios. This is the subject of 
this paper. 

Throughout we have taken go = and Hq = 
100 h km s^^Mpc"^ 



2 METHOD 

2.1 The relationship between galaxy and AGN 
luminosity functions 

The luminosity function of AGN must be related to the 
mass function of massive black holes. If the masses of cen- 
tral black holes are related to the masses of their host 
galaxy bulges (Magorrian et al. 1998, Kormendy & Rich- 
stone 1995), then the luminosity functions of AGN and 
galaxies must be strongly related. This is modeled in the 
following formalism. 

We define the luminosity function as 

~ dVdL 

Where is number of objects, V is comoving volume and 
L is luminosity, and start from the luminosity functions ipi 
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Table 1. Schochtor model galaxy luminosity functions from 
Folkos et al. (1999) and assumed fraction fsph of B band light 
due to the spheroidal component. 



Type 


M* 


a 




fsph 


E/SO 


-19.61 


-0.74 


9.0 


0.7 


Sab 


-19.68 


-0.86 


3.9 


0.25 


Sbc 


-19.38 


-0.99 


5.3 


0.15 


Scd 


-19.00 


-1.21 


6.5 


0.1 


Sdm 


-19.02 


-1.73 


2.1 


0.02 



of galaxies of n different morphological types i. We assume 
that the spheroid components of galaxies of type i produce 
a fraction of their light fsph{i)- We can then translate from 
the (j>iS to a spheroid luminosity function (j)sph by: 

(j>i{Lsph/ fsphii)) 



4'sph{Lsph) 



1=1 



fsph{i) 



We then obtain the mass function of spheroids by: 



dVdMsph dMsph 

and a black hole mass function <^bh is obtained using 
d'^N f d'^N 

where P{Mbh \ M^ph) is the distribution of nuclear black 
hole masses given spheroid mass Msph- 

Some fraction /agn of nuclear black holes are in 
an active (luminous) state and emit with luminosity 
Lagn- Both Lagn and Jagn are likely to depend on 
Mbh- For simplicity we assume a functional relationship 
Lagn ~ F[Mbh), although a distribution of luminosities 
P{Lagn I Mbh) would be more realistic. 

The AGN luminosity function <j)AGN is then: 

dMBH 



4>AGN = JaGN^BH- 



(1) 



dLAGN 

It is appropriate to formulate the ratio of Lagn 
to Mbh in terms of the Eddington ratio e, the ra- 
tio of bolometric to Eddington luminosity Le = 1.3 x 
10^^ {Mbh /MQ)eTgs-\ 
_ Lagn Mq 



Mbh 1.3 x lO^Serg s' 



(2) 



3 CONSTRUCTING THE BLACK HOLE 
MASS FUNCTION 

For the local galaxy luminosity functions we have chosen 
to use the recent determinations from the 2DF galaxy red- 
shift survey (Folkes et al. 1999). This provides Schechter 
function model luminosity functions for 5 different spec- 
troscopic classes of galaxy (corresponding to different mor- 
phological types) at 2 < 0.2. We have used values for fsph 
(the fraction of B band light due to the spheroid com- 
ponent) based on the values given by Meisels & Ostriker 
(1984). 

The model luminosity functions and spheroid frac- 
tions are listed in Table 0. 

For the spheroid mass - to light ratio we use the best 
fit relation of Magorrian et al. (1998) for V: 
M/Mq = 0.097 h (L/Lq)'--''^ 

and assume B-V=0.9 for the spheroid component of all 
galaxy types (Pence 1976). For the black hole to spheroid 




log M (h \) 
Figure 1. Mass function of black holes derived in Section 



mass distribution we use the best fit log-gaussian relation 
from Magorrian et al. (1998): 

PiMBH I Msph) = — 

MsifA^loge(lO) 

where A = 0.51, /bh = 0.97 and xo = 5.2 x 10"^ We 
truncate this distribution at ±3(t because otherwise it pre- 
dicts an unrealistically large population of very massive 
(> IO^^Mq) black holes. 

The resultant black hole mass function is shown in 
Fig. §. 



4 FITTING THE ACCRETION RATE AND 
ACTIVE FRACTION 

For the AGN luminosity function we have chosen to use 
the extended Etnstem medium sensitivity survey (EMSS) 
X-ray selected sample (Stocke et al. 1991) restricted to 
type 1 (broad line) Seyferts with z < 0.2. X-ray selection 
is ideal because it is insensitive to the optical host galaxy 
properties. The restriction to broad line objects ensures 
that the sample contains only unobscured AGN for which 
the observed X-ray flux is a good measure of intrinsic 
luminosity, and the redshift restriction ensures that the 
sample is well matched to the 2DF galaxy sample and is 
without significant cosmological evolution. 

The simplest relation between black hole mass and 
AGN luminosity (and therefore the starting point for this 
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Table 2. The nine closest galaxies from the EMSS Seyfert 1 
sample, with their de Vaucouleurs (1959) morphological classi- 
fications as listed in the NASA Extragalactic Database (NED). 



Source 


logLx 


morpholog; 


MS1215.9+3005 


2.18 


SB(s)a 


MS0339.8-2124 


2.30 


SB(rs)a 


MS0459. 5+0327 


2.59 


E 


MS1158. 6-0323 


2.45 


E 


MS2252. 2-1-1126 


2.00 


Sab 


MS1846.5-7857 


2.08 


SAB(r)b 


MS1136.5-I-3413 


2.74 


SBO 


MS0048.8-I-2907 


2.96 


SB(s)b 



investigation) is obtained by assuming an Eddington ra- 
tio e and an active fraction /agn independent of mass or 
galaxy type; this will be referred to as Model A henceforth. 
Since we are constructing an X-ray luminosity function 
rather than a bolometric luminosity function, and only 
considering type 1 Seyferts rather than all AGN, we re- 
place fAGN with fsi and reformulate Equation ^ in terms 
of the Eddington ratio and include a bolometric correction 
/3, where /3 is the ratio of the bolometric luminosity to the 
X-ray luminosity. 
O.OlSeMBH 

= -j^ (3) 

where Lx is the AGN 0.3-3.5 keV luminosity in units of 
10''°ergs~^ and Mbh is in units of Mq. 

'''^ dVdLx O.OlSe dVdMsH 

We adopt /3 = 20 based on the mean type 1 AGN spectral 
energy distribution from Elvis et al. (1994). 

A binned (AlogL = 0.3) AGN X-ray luminosity 
function was computed from the EMSS sample (Fig. 2a) 
and compared to the model predicted luminosity function 
using the method described in Sections 2.3 and 5.2 of Page 
& Carrera (2000). Only luminosity bins containing > 10 
objects were used in the fitting so that could be used 
as a goodness of fit estimator; the resultant binned lumi- 
nosity function uses 152 EMSS AGN. 

A very good fit to the X-ray luminosity function of 
type 1 AGN is easily found. The best fit has he — 1.8 ± 
0.5 X 10"^ and fsi = 7 ± 2 x 10"^ (where the errors 
define a box containing the la, Ax^ = 2.3, confidence 
interval) with an extremely low x^ of 0.4 for 6 fitted data 
points and 2 free parameters (i.e. 4 degrees of freedom, so 
X^//^ ~ 0.1). The best fit luminosity function is shown in 
Fig. ^; the predicted luminosity function passes almost 
exactly through the data. Fig. Hb shows Ax^ confidence 
contours on the fitted parameters. 

However, despite the very good fit, there is a prob- 
lem with this simple model. It is seen in Fig. ^ that the 
EMSS luminosity function is produced almost exclusively 
by E and SO galaxies with only a small contribution from 
Sab galaxies at the low luminosity end, but it is well known 
that Seyferts are often spirals. Indeed, spirals are a signif- 
icant component of the EMSS Seyfert sample itself; this 
is demonstrated in Table ^ which gives the morphologies 
of the nearest EMSS Seyferts. 

To reproduce the the EMSS luminosity function with 
a mix of morphological types requires additional complex- 
ity in the model: to prevent early type galaxies completely 
dominating the luminosity function, fsi and/or e must de- 



pend on Mbh and/or galaxy morphology. Since different 
galaxy morphologies dominate the mass function at differ- 
ent masses, a direct dependency of fsi or e on Mbh results 
in an indirect dependency on morphology and vice-versa. 
If e is to be varied it must be in such a way as to steepen 
the model luminosity function so that the relative contri- 
bution of E/SO galaxies is reduced in the EMSS luminosity 
range. This means that more massive black holes (in E/SO 
galaxies) must accrete with lower Eddington ratios than 
those of lower mass. If fsi is to be varied it must be such 
that a smaller proportion of more massive black holes are 
actively accreting Seyferts than those of lower mass. 

An ad-hoc example of one of these models is shown 
in Figs. In this model (hereafter Model B) we have as- 
sumed that the fraction of active Seyferts in E/SO galax- 
ies, and their Eddington ratios, are only half that of later 
type galaxies. Again, the best fit luminosity function is 
an extremely good fit with a x^ oi 0.45, but this time 
the AGN luminosity function is produced by AGN with a 
more plausible mix of host galaxy morphologies. The fitted 
parameters are fsi = 7 ± 3 x 10~^ and he — ^ ±2 x 10"^ 
for Sab and later galaxies (and by design half these values 
for E and SO galaxies). 

To obtain a conservative limit on the range of param- 
eter space that brackets fsi and e, we have also considered 
the opposite extreme to the E/SO dominated Model A: in 
Model C the active fraction of E and SO galaxies is set 
to zero and they therefore make no contribution whatso- 
ever to the AGN luminosity function. Since the fraction 
of Seyfert nuclei which are hosted by E and SO galaxies is 
certainly between and 1, it is reasonable to expect that 
the real values of fsi and e must lie somewhere between 
the acceptable values for Model A and Model C. The best 
fit AGN luminosity function and Ax^ confidence contours 
for Model C are shown in Fig. 0. 



5 DISCUSSION 

The approach used here easily reproduces the shape of 
the AGN luminosity function; models A, B and C all have 
extremely low x^- The steepening of the AGN luminosity 
function from low to high luminosity is inherited from the 
black hole mass function, which in turn inherits the shape 
from the galaxy luminosity functions. Hence the shape of 
the AGN luminosity function ultimately derives from the 
same physical processes that give rise to the shapes of 
galaxy luminosity functions. 

The fitted Seyfert 1 fraction in Sab and later type 
galaxies for Model B is ~ 0.5 — 1%. This is quite a robust 
result: it isn't strongly affected by the contribution of E 
and SO galaxies because Models A and C have very similar 
best fit values for fsi - This compares well with the findings 
of optical emission line surveys of local galaxies. Maiolino 

6 Rieke (1995), find that 5% of the revised Shapley-Ames 
catalogue of galaxies are Seyferts (with a Seyfert 1: Seyfert 
2 ratio of 1:4 this corresponds to a Seyfert 1 fraction of 
1%). Ho, Filippenko & Sargent (1997), find a Seyfert frac- 
tion which is twice that found by Maiolino & Rieke (1995), 
and therefore (assuming the Maiolino & Rieke Seyfert 1: 
Seyfert 2 ratio) about a factor two higher than would be 
expected from our value of fsi- However, their Seyfert 
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log hyr (h ^ 10*° erg s ^ 

Figure 3. (a) X— ray luminosity function of the EMSS z < 0.2 AGN and best fit model in which the Seyfert 1 fraction of E and SO 
galaxies is only half that in other types, and the Seyfert 1 nuclei hosted by E and SO galaxies accrete with only half the Eddington 
ratio of Seyfert nuclei in later galaxy types (Model B). (b) confidence contours for the values of Eddington ratio e and Seyfert 1 
fraction fgi for Sab and later type galaxies in Model B. 



detection rate is high because their survey is sensitive to 
objects with very weak emission lines, which they term 
"dwarf" Seyferts; it is not yet known whether these ob- 
jects have the same X-ray properties as more luminous 
Seyferts which make up the EMSS sample. Both Maiolino 
& Rieke (1995) and Ho, Filippenko & Sargent (1997) find 
that the most common morphology for Seyfert galaxies 
is Sa - Sb further justifying our rejection of model A in 
Section 13. 



In contrast to /si, the value of e for typical Sab 
Seyfert galaxies depends quite strongly on the Seyfert ac- 
tivity in E and SO galaxies. Taking the extremes of mod- 
els A and C, a conservative conclusion is that 0.001 < 
e < 0.02 for Seyfert nuclei in Sab galaxies. Our reasonable 
model B suggests e ~ 0.005 should be typical, consider- 
ably sub-Eddington and consistent with Seyfert activity 
occurring recurrently in a significant fraction of galaxies 
(model 'R' from Cavaliere & Padovani 1988). This is in 
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Figure 4. (a) X— ray luminosity function of the EMSS z < 0.2 AGN and best fit model in which E and SO galaxies host no AGN 
(Model C). (b) confidence contours for the values of Eddington ratio e and Seyfert 1 fraction /gi for Sab and later type galaxies 
in Model C. 



good agreement with the predictions of recent models in 
which AGN, with lifetimes of order 10^ years, are pro- 
duced during the formation and merging of galaxies in 
a cold dark matter dominated Universe; Kauffmann & 
Haehnelt (2000) predict typical e ~ 0.01 while Haiman 
& Menou (2000) predict typical e ~ 0.001. Our results are 
also consistent with the range of Seyfert Eddington ratios 
given in the recent compilation by Wandel (1999). These 
lie between 0.001 and 1 and were obtained using a variety 
of different methods: broad line region kinematics. X-ray 
variability and modelling accretion disk spectra. 

However, the range of acceptable e is too low to be 
consistent with accretion disk models which have large 
outbursts. Siemiginowska & Elvis (1997) show that if AGN 
accretion disks are subject to the thermal-viscous insta- 
bility driven by hydrogen ionization (Lin & Shields 1986) 
they are probably only observed in outburst, with typi- 
cal e ~ 0.1. The study by Burderi, King & Szuszkiewicz 
(1998) concludes that if AGN have optically thick, geo- 
metrically thin accretion disks (Shakura & Sunyaev 1973) 
they almost certainly are subject to this kind of insta- 
bility. The results presented here would therefore suggest 
that AGN are not fuelled by standard thin disks. 

As explained in Section ^, the only way to get a 
plausible distribution of Seyfert host galaxies is to model 
the more massive black holes in early type galaxies with 
smaller values of e, fsi or both. Both these alternatives 
have important consequences for our understanding of 
AGN behaviour. 

The first (e smaller for more massive AGN) is intu- 
itively reasonable because the early type galaxies which 
host more massive AGN contain less gas with which to 
feed them. However, it is contrary to the results of Wan- 
del & Petrosian (1988) and Sun & Malkan (1989), who 
find from accretion disk modelling that luminous QSOs 
are both more massive, and accrete with higher Edding- 



ton ratios, than Seyfert galaxies (although this could be 
interpreted as a dependence of Eddington ratio on red- 
shift rather than mass). The AGN evolution models of 
Kauffmann & Haehnelt (2000) and Salucci et al. (1999), 
also favour a situation in which the more luminous present 
epoch AGN are more massive and have higher Eddington 
ratios. 

The second way to obtain reasonable Seyfert mor- 
phologies, by modelling a smaller Seyfert fraction fsi for 
more massive black holes, is less controversial. For exam- 
ple, both Maiolino & Rieke (1995) and Ho, Filippenko & 
Sargent (1997) find that a higher fraction of Sa-Sb galax- 
ies have Seyfert nuclei than E/SO galaxies. In this case 
the mass function of active nuclei must be steeper than 
the mass function of inactive nuclei, unlike the simple re- 
current activity models represented by model 'R' in figure 
1 of Cavahere & Padovani (1988). 

This investigation can be substantially improved 
upon with a sample of X-ray selected Seyferts with opti- 
cal morphologies; this would allow the matching of galaxy 
and AGN luminosity functions individually for each mor- 
phological type, providing a much more rigorous solution 
than our model B. This obviously requires a larger AGN 
sample than the EMSS one used here, for which the lumi- 
nosity functions of different morphological types would be 
dominated by Poisson noise. We can expect that such a 
sample may soon be available from ROSAT All Sky Survey 
optical identification programmes. 
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